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Abstract. The supermassive black holes observed at the centers of almost all present- 
day galaxies, had a profound impact on their environment. I highlight the principle of 
self-regulation, by which supermassive black holes grow until they release sufficient en- 
ergy to unbind the gas that feeds them from their host galaxy. This principle explains 
several observed facts, including the correlation between the mass of a central black hole 
and the depth of the gravitational potential well of its host galaxy, and the abundance 
and clustering properties of bright quasars in the redshift interval of z ~ 2-6. At lower 
redshifts, quasars might have limited the maximum mass of galaxies through the sup- 
pression of cooling flows in X-ray clusters. The seeds of supermassive black holes were 
likely planted in dwarf galaxies at redshifts z > 10, through the collapse of massive or 
supermassive stars. The minimum seed mass can be identified observationally through 
the detection of gravitational waves from black hole binaries by Advanced LIGO or 
LISA. Aside from shaping their host galaxies, quasar outflows filled the intergalactic 
medium with magnetic fields and heavy elements. Beyond the reach of these outflows, 
the brightest quasars at z > 6 have ionized exceedingly large volumes of gas (tens of 
comoving Mpc) prior to global reionization, and must have suppressed the faint end of 
the galaxy luminosity function in these volumes before the same occurred through the 
rest of the universe. 



1 The Principle of Self-Regulation 

The fossil record in the present-day universe indicates that every bulged galaxy 
hosts a supermassive black hole (BH) at its center pQ. These BHs are dormant 
or faint most of the time, but ocassionally flash in a short burst of radiation that 
lasts for a small fraction of the Hubble time. The short duty cycle acounts for 
the fact that bright quasars are much less abundant than their host galaxies, but 
it begs the more fundamental question: why is the quasar activity so brief? A 
natural explanation is that quasars are suicidal, namely the energy output from 
the BHs regulates their own growth. 

Supermassive BHs make up a small fraction, < 10 -3 , of the total mass in their 
host galaxies, and so their direct dynamical impact is limited to the central star 
distribution where their gravitational influence dominates. Dynamical friction on 
the background stars keeps the BH close to the center. Random fluctuations in 
the distribution of stars induces a Brownian motion of the BH. This motion can 
be decribed by the same Langevin equation that captures the motion of a massive 
dust particle as it responds to random kicks from the much lighter molecules of 
air around it PJ. The characteristic speed by which the BH wanders around the 
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center is small, ~ (m^/MBu) 1 ^ 2 ^*, where m* and -Mbh are the masses of a single 
star and the BH, respectively, and cr* is the stellar velocity dispersion. Since the 
random force fluctuates on a dynamical time, the BH wanders across a region 
that is smaller by a factor of ~ {iti^/Mbb) 1 ^ 2 than the region traversed by the 
stars inducing the fluctuating force on it. 

The dynamical insignificance of the BH on the global galactic scale is mis- 
leading. The gravitational binding energy per rest-mass energy of galaxies is 
of order ~ (<t*/c) 2 < 10 -6 . Since BH are relativistic objects, the gravitational 
binding energy of material that feeds them amounts to a substantial fraction its 
rest mass energy. Even if the BH mass occupies a fraction as small as ~ 10 -4 
of the baryonic mass in a galaxy, and only a percent of the accreted rest-mass 
energy leaks into the gaseous environment of the BH, this slight leakage can 
unbind the entire gas reservoir of the host galaxy! This order-of-magnitude es- 
timate explains why quasars are short lived. As soon as the central BH accretes 
large quantities of gas so as to significantly increase its mass, it releases large 
amounts of energy that would suppress further accretion onto it. In short, the 
BH growth is self-regulated. 

The principle of self-regulation naturally leads to a correlation between the 
final BH mass, Mbh, and the depth of the gravitational potential well to which 
the surrounding gas is confined, ~ cr 2 . Indeed such a correlation is observed in 
the present-day universe |2j • The observed power-law relation between ATbh and 
cr* can be generalized to a correlation between the BH mass and the circular 
velocity of the host halo, v c which in turn can be related to the halo mass, 
■Mhaio, and redshift, z ^B] 

Af bh (A/haio, z) = const x v 5 c 

= e M hal0 ( I ^^) 3 [C(z)]l(l + ,)i, (1) 

where e « 10 -5 7 is a constant, ((z) is close to unity and defined as ( = 
[(n m /Q^)(A c /18n% f2* m = [l + ([2 A /f2 m )(l + z)- 3 }-\ A c = 187r 2 + 82d-39d 2 , 
and d = — 1 (see equations 22-25 in Ref. [H] for the relation between v c and 
■Mhalo)- If quasars shine near their Eddington limit as suggested by observations 
of low and high-redshift quasars jril7j . then the above value of e implies that 
a fraction of ~ 5-10% of the energy released by the quasar over a galactic 
dynamical time needs to be captured in the surrounding galactic gas in order 
for the BH growth to be self-regulated JH] . 

With this interpretation, the Mbb-0* relation reflects the limit introduced 
to the BH mass by self-regulation; deviations from this relation are inevitable 
during episodes of BH growth or as a result of mergers of galaxies that have no 
cold gas in them. A physical scatter around this upper envelope could also result 
from variations in the efficiency by which the released BH energy couples to the 
surrounding gas. 

Various prescriptions for self-regulation were sketched by Silk & Rees [H] . 
These involve either energy or momentum-driven winds, where the latter type 
is a factor of ~ v c /c less efficient [13114115] . Wyithe & Loeb JHI demonstrated 
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Fig. 1. Comparison of the observed and model luminosity functions (from The 
data points at z < 4 are summarized in Ref. 0, while the light lines show the double 
power-law fit to the 2dF quasar luminosity function At z ~ 4.3 and z ~ 6.0 
the data is from Refs. |10) . The grey regions show the 1-cr range of logarithmic slope 
([—2.25, —3.75] at 2 ~ 4.3 and [—1.6, —3.1] at z ~ 6), and the vertical bars show the 
uncertainty in the normalization. The open circles show data points converted from 
the X-ray luminosity function |llj of low luminosity quasars using the median quasar 
spectral energy distribution. In each panel the vertical dashed lines correspond to the 
Eddington luminosities of BHs bracketing the observed range of the Mbh~«c relation, 
and the vertical dotted line corresponds to a BH in a 10 13 5 M© galaxy. 



that a particularly simple prescription for an energy-driven wind can reproduce 
the luminosity function of quasars out to highest measured redshift, z ~ 6 (see 
Figs- H] and as well as the observed clustering properties of quasars at z ~ 3 
[T7] (see Fig.^J). The prescription postulates that: (i) self- regulation leads to the 
growth of Mbh up the redshift-indcpendent limit as a function of v c in Eq. Q , 
for all galaxies throughout their evolution; and (ii) the growth of Mbh to the 
limiting mass in Eq. occurs through halo merger episodes during which the 
BH shines at its Eddington luminosity (with the median quasar spectrum) over 
the dynamical time of its host galaxy, idyn- This model has only one adjustable 
parameter, namely the fraction of the released quasar energy that couples to the 
surrounding gas in the host galaxy. This parameter can be fixed based on the 
Mbh-o* relation in the local universe 0]. It is remarkable that the combination 
of the above simple prescription and the standard /1CDM cosmology for the 
evolution and merger rate of galaxy halos, lead to a satisfactory agreement with 
the rich data set on quasar evolution over cosmic history. 
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The cooling time of the heated gas is typically longer than its dynamical time 
and so the gas should expand into the galactic halo and escape the galaxy if its 
initial temperature exceeds the virial temperature of the galaxy |lfi| . The quasar 
remains active during the dynamical time of the initial gas reservoir, ~ 10 7 years, 
and fades afterwards due to the dilution of this reservoir. Accretion is halted as 
soon as the quasar supplies the galactic gas with more than its binding energy. 
The BH growth may resume if the cold gas reservoir is replenished through a 
new merger. 




Fig. 2. The comoving density of supermassive BHs per unit BH mass (from |16|1. The 
grey region shows the estimate based on the observed velocity distribution function of 
galaxies in Ref. |18| and the M^-Vc relation in Eq. (0. The lower bound corresponds 
to the lower limit in density for the observed velocity function while the grey lines 
show the extrapolation to lower densities. We also show the mass function computed 
at z = 1, 3 and 6 from the Press-SchechterfH?) halo mass function and Eq. Q, as 
well as the mass function at z ~ 2.35 and z ~ 3 implied by the observed density of 
quasars and a quasar lifetime of order the dynamical time of the host galactic disk, 
tdyn (dot-dashed lines). 

Agreement between the predicted and observed correlation function of quasars 
(Fig. I3J) is obtained only if the BH mass scales with redshift as in Eq. (JTJ and 
the quasar lifetime is of the order of the dynamical time of the host galactic disk 

tdyn = io 7 [£(z)]- 1/2 (^) yr. (2) 

The inflow of cold gas towards galaxy centers during the growth phase of the 
BH would naturally be accompanied by a burst of star formation. The fraction 
of gas that is not consumed by stars or ejected by supernovae, will continue to 
feed the BH. It is therefore not surprising that quasar and starburst activities 
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co-exist in Ultra Luminous Infrared Galaxies |21| . and that all quasars show 
broad metal lines indicating a super-solar metallicity of the surrounding gas |22| . 
Applying a similar self-regulation principle to the stars, leads to the expectation 
|lfil2.'-i| that the ratio between the mass of the BH and the mass in stars is 
independent of halo mass (as observed locally [21]) but increases with redshift 
as oc ^(z) 1 ^ 2 (l + z) 3 / 2 . A consistent trend has indeed been inferred in an observed 
sample of gravitationally-lensed quasars [25] . 




Fig. 3. Predicted correlation function of quasars at various redshifts in comparison to 
the 2dF data |2U| (from |17|1. The dark lines show the correlation function predictions 
for quasars of various apparent B-band magnitudes. The 2dF limit is B ~ 20.85. The 
lower right panel shows data from entire 2dF sample in comparison to the theoretical 
prediction at the mean quasar redshift of (z) = 1.5. The B = 20.85 prediction at this 
redshift is also shown by thick gray lines in the other panels to guide the eye. The 
predictions are based on the scaling A/bh oc Vc in Eq. Q . 
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The upper mass of galaxies may also be regulated by the energy output 
from quasar activity. This would account for the fact that cooling flows are 
suppressed in present-day X-ray clusters |2fil27l28| . and that massive BHs and 
stars in galactic bulges were already formed at z ~ 2. The quasars discovered 
by the Sloan Digital Sky Survey (SDSS) at z ~ 6 mark the early growth of 
the most massive BHs and galactic spheroids. The present-day abundance of 
galaxies capable of hosting BHs of mass ~ 10 9 M Q (based on Eq. 0) already 
existed at z ~ 6 |29| . At some epoch, the quasar energy output may have led to 
the extinction of cold gas in these galaxies and the suppression of further star 
formation in them, leading to an apparent "anti-hierarchical" mode of galaxy 
formation where massive spheroids formed early and did not make new stars 
at late times. In the course of subsequent merger events, the cores of the most 
massive spheroids acquired an envelope of collisionless matter in the form of 
already- formed stars or dark matter [29| . without the proportional accretion of 
cold gas into the central BH. The upper limit on the mass of the central BH and 
the mass of the spheroid is caused by the lack of cold gas and cooling flows in 
their X-ray halos. In the cores of cooling X-ray clusters, there is often an active 
central BH that supplies sufficient energy to compensate for the cooling of the 
gas |27l2fill3| . The primary physical process by which this energy couples to the 
gas is still unknown. 

2 Feedback on Large Intergalactic Scales 

Aside from affecting their host galaxy, quasars disturb their large-scale cosmolog- 
ical environment. Powerful quasar outflows are observed in the form of radio jets 
|3()| or broad-absorption-line winds j31| . The amount of energy carried by these 
outflows is largely unknown, but could be comparable to the radiative output 
from the same quasars. Furlanetto & Loeb [32] have calculated the intergalactic 
volume filled by such outflows as a function of cosmic time (see Fig. QJ. This 
volume is likely to contain magnetic fields and metals, providing a natural source 
for the observed magnetization of the metal-rich gas in X-ray clusters [S3J and 
in galaxies |34| . The injection of energy by quasar outflows may also explain the 
deficit of Lya absorption in the vicinity of Lyman-break galaxies 35 36^ and the 
required pre- heating in X-ray clusters |37I27| . 

Beyond the reach of their outflows, the brightest SDSS quasars at z > 6 are 
inferred to have ionized exceedingly large regions of gas (tens of comoving Mpc) 
around them prior to global reionization (see Fig. and Refs. 38 430- Thus, 
quasars must have suppressed the faint-end of the galaxy luminosity function in 
these regions before the same occurred throughout the universe. The recombi- 
nation time is comparable to the Hubble time for the mean gas density at z ~ 7 
and so ionized regions persist |39j on these large scales where inhomogeneities 
are small. The minimum galaxy mass is increased by at least an order of mag- 
nitude to a virial temperature of ~ 10 5 K in these ionized regions It would 
be particularly interesting to examine whether the faint end (cr* < 30km s _1 ) of 




Fig. 4. The global influence of magnetized quasar outflows on the intergalactic medium 
(from 32 ). Upper Panel: Predicted volume filling fraction of magnetized quasar bub- 
bles F(z), as a function of redshift. Lower Panel: Ratio of normalized magnetic energy 
density, ug/e-i, to the fiducial thermal energy density of the intergalactic medium 
Ufid = 3n(z)kTiGM , where Tigm = 10 4 K, as a function of redshift (see \62\ for 
more details). In each panel, the solid curves assume that the blast wave created by 
quasar ouflows is nearly (80%) adiabatic, and that the minimum halo mass of galax- 
ies, Mh, m in, is determined by atomic cooling before reionization and by suppression 
due to galactic infall afterwards (top curve), Mh : min = W 9 Mq (middle curve), and 

min — 

1O 1O M (bottom curve). The dashed curve assumes a fully-radiative blast 
wave and fixes Mh,min by the thresholds for atomic cooling and infall suppression. The 
vertical dotted line indicates the assumed redshift of complete reionization, z r = 7. 



the luminosity function of dwarf galaxies shows any moduluation on large-scales 
around rare massive BHs, such as M87. 

To find the volume filling fraction of relic regions from z ~ 6, we consider 
a BH of mass Mbh ~ 3 x 1O 9 M . We can estimate the co- moving density of 



8 Avi Loeb 



BHs directly from the observed quasar luminosity function and our estimate of 
quasar lifetime. At z ~ 6, quasars powered by Mbh ~3x 10 9 M Q BHs had a 
comoving density of ~ 0.5Gpc _3 5Ej- However, the Hubble time exceeds idyn 
by a factor of ~ 2 x 10 2 (reflecting the square root of the overdensity in cores 
of galaxies), so that the comoving density of the bubbles created by the z ~ 6 
BHs is ~ 10 2 Gpc~ 3 (see Fig. |5J). The density implies that the volume filling 
fraction of relic z ~ 6 regions is small, < 10%, and that the nearest BH that had 
Mbh ~3x 10 9 M Q at z ~ 6 (and could have been detected as an SDSS quasar 

1/3 

then) should be at a distance dbh ~ (47r/3 x 10 2 ) Gpc ~ 140Mpc which is 
almost an order-of-magnitude larger than the distance of M87, a galaxy known 
to possess a BH of this mass |4"3) . 

What is the most massive BH that can be detected dynamically in a local 
galaxy redshift survey? SDSS probes a volume of ~ lGpc 3 out to a distance 
~ 30 times that of M87. At the peak of quasar activity at z ~ 3, the density 
of the brightest quasars implies that there should be ~ 100 BHs with masses of 
3 x 1O 1O M per Gpc 3 , the nearest of which will be at a distance c?bh ~ 130Mpc, 
or ~ 7 times the distance to M87. The radius of gravitational influence of the BH 
scales as M hh /v% oc M^ 5 . We find that for the nearest 3 x 1O 9 M and 3 x 10 10 M Q 
BHs, the angular radius of influence should be similar. Thus, the dynamical 
signature of ~ 3 x 10 M BHs on their stellar host should be detectable. 



3 What seeded the growth of the supermassive black 
holes? 

The BHs powering the bright SDSS quasars possess a mass of a few x 1O 9 M , and 
reside in galaxies with a velocity dispersion of ~ 500km s _1 45 . A quasar radiat- 
ing at its Eddington limiting luminosity, L E = 1.4 x 10 46 erg s _1 (M b h/TO 8 M ), 
with a radiative efficiency, e ra d = Le/Mc 2 would grow exponentially in mass 
as a function of time t, Mbh = M see d exp{t/ie} on a time scale, ts = 4.1 x 
10 7 yr(e ra d/0.1). Thus, the required growth time in units of the Hubble time 
thubbie = 9 x 10 8 yr[(l + z)/7]- 3 / 2 is 

growth / £rad \( 1 + Z V'* ( M hb /10 9 M Q \ 

^hubblc 

uo%7 V 7 / V A/ sood/iooM y ■ w 

The age of the universe at z ~ 6 provides just sufficient time to grow an SDSS 
BH with Mbh ~ 1O 9 M out of a stellar mass seed with e ra d = 10% The 
growth time is shorter for smaller radiative efficiencies, as expected if the seed 
originates from the optically-thick collapse of a supermassive star (in which case 
M see d in the logarithmic factor is also larger). 

What was the mass of the initial BH seeds? Were they planted in early dwarf 
galaxies through the collapse of massive, metal free (Pop-Ill) stars (leading to 
M see d of hundreds of solar masses) or through the collapse of even more mas- 
sive, i.e. supermassive, stars \47Jj ? Bromm & Loeb [48] have shown through 
a hydrodynamical simulation (see Fig. [SJ that supermassive stars were likely to 
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Fig. 5. Quasars serve as probes of the end of reionization. The measured size of the 
HII regions around SDSS quasars can be used |4(JI41| to demonstrate that a significant 
fraction of the intergalactic hydrogen was neutral at z ~ 6.3 or else the inferred size of 
the quasar HII regions would have been much larger than observed (assuming typical 
quasar lifetimes |42|V Also, quasars can be used to measure the redshift at which 
the intergalactic medium started to transmit Lya photonspfe 43 . The upper panel 
illustrates how the line-of-sight towards a quasar intersects this transition redshift. The 
resulting Lya transmission of the intrinsic quasar spectrum is shown schematically in 
the lower panel. 



form in early galaxies at z ~ 10 in which the virial temperature was close to the 
cooling threshold of atomic hydrogen, ~ 10 4 K. The gas in these galaxies con- 
densed into massive ~ 10 6 M Q clumps (the progenitors of supermassive stars), 
rather than fragmenting into many small clumps (the progenitors of stars), as 
it does in environments that are much hotter than the cooling threshold. This 
formation channel requires that a galaxy be close to its cooling threshold and 
immersed in a UV background that dissociates molecular hydrogen in it. These 
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Fig. 6. SPH simulation of the collapse of an early dwarf galaxy with a virial temper- 
ature just above the cooling threshold of atomic hydrogen and no H2 (from |48|1. The 
image shows a snapshot of the gas density distribution at z ~ 10, indicating the forma- 
tion of two compact objects near the center of the galaxy with masses of 2.2 x 1O 6 M0 
and 3.1 x 10 6 Mq, respectively, and radii < 1 pc. Sub-fragmentation into lower mass 
clumps is inhibited as long as molecular hydrogen is dissociated by a background UV 
flux. These circumstances lead to the formation of supermassive stars |47| that in- 
evitably collapse and trigger the birth of supermassive black holes 47 49 . The box size 
is 200 pc. 

requirements should make this channel sufficiently rare, so as not to overproduce 
the cosmic mass density of supermassive BH. 

The minimum seed BH mass can be identified observationally through the 
detection of gravitational waves from BH binaries with Advanced LIGO [5TJ or 
with LISA [5U| . Most of the mHz binary coalescence events originate at z > 7 
if the earliest galaxies included BHs that obey the M^-Vc relation in Eq. (JIJ. 
The number of LISA sources per unit redshift per year should drop substantially 
after reionization, when the minimum mass of galaxies increased due to photo- 
ionization heating of the intergalactic medium. Studies of the highest redshift 
sources among the few hundred detectable events per year, will provide unique 
information about the physics and history of BH growth in galaxies ■ 

The early BH progenitors can also be detected as unresolved point sources, 
using the future James Webb Space Telescope (JWST). Unfortunately, the spec- 
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trum of metal-free massive and supermassive stars is the same, since their surface 
temperature ~ 10 5 K is independent of mass Hence, an unresolved cluster 
of massive early stars would show the same spectrum as a supermassive star of 
the same total mass. 

In closing, it is difficult to ignore the possible environmental impact of quasars 
on anthropic selection. One may wonder whether it is not a coincidence that our 
Milky- Way Galaxy has a relatively modest BH mass of only a few million solar 
masses in that the energy output from a much more massive (e.g. ~ 10 9 A/q) 
black hole would have disrupted the evolution of life on our planet. A proper 
calculation remains to be done (as in the context of nearby Gamma-Ray Bursts 
|54|) in order to demonstrate any such link. 
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